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Summary 

This  paper  reconsiders  the  problem  of  determining  the  elastostatic 
field  near  the  tip  of  a  crack  in  an  all-around  infinite  body  deformed  by  a 
"Mode  III"  loading  at  infinity  to  a  state  of  anti-plane  shear.  The  problem 
is  treated  for  a  class  of  incompressible,  homogeneous,  isotropic  elastic 
materials  whose  constitutive  laws  permit  a  loss  of  ellipticity  in  the  govern¬ 
ing  displacement  equation  of  equilibrium  at  sufficiently  severe  shearing 
strains.  The  analysis  represents  a  generalization  of  that  reported  in  an 
earlier  study  and,  as  before,  is  carried  out  for  the  "small-scale  nonlinear 
crack  problem",  in  which  a  crack  of  finite  length  is  replaced  by  a  semi¬ 
infinite  one,  and  the  nonlinear  field  far  from  the  crack-tip  is  matched  to 
the  near  field  predicted  by  the  linearized  theory.  The  methods  employed 
in  the  present  paper  are  necessarily  largely  qualitative,  since  they  apply 
to  all  materials  in  the  class  considered.  The  principal  feature  of  the 
resulting  elastic  field  is  the  presence  of  two  symmetrically  located  curves 


issuing  from  the  crack-tip  and  bearing  discontinuities  in  displacement 


The  results  communicated  in  this  paper  were  obtained  in  the  course  of  an 
investigation  supported  in  part  by  Contract  N00014-75-C-0196  with  the 
Office  of  Naval  Research  in  Washington,  D.  C. 


gradient  and  stress.  \ 

K 

Introduction 

In  a  recent  paper  [  l]  we  investigated  the  elastostatic  field  near 
the  tip  of  a  crack  under  conditions  of  finite  anti-plane  shear  for  a  very 
special  hypothetical  incompressible  —  homogeneous  and  isotropic  —  elastic 
material  whose  constitutive  law  permits  a  loss  of  ellipticity  of  the  appropriate 
displacement  equation  of  equilibrium  at  sufficiently  severe  deformations. 

The  problem  treated  in  [  l]  involves  an  all-around  infinite  body  containing 
a  plane  crack  of  constant  width  and  infinite  length.  At  infinity  the  body  is 
subjected  to  a  state  of  simple  shear  parallel  to  the  edges  of  the  crack.  The 
analysis  reported  in  [  l]  is  carried  out  within  the  framework  of  finite  elas¬ 
ticity  theory. 

When  the  amount  of  shear  applied  at  infinity  is  small,  the  field  near 
either  crack-tip  can  be  determined  on  the  basis  of  an  asymptotic  scheme 
in  which  the  crack  of  finite  width  is  replaced  by  one  of  semi -infinite  width, 
and  the  far  field  is  required  to  match  the  elastostatic  field  near  the  crack- 
tip  predicted  by  the  solution  to  the  original  problem  according  to  the  linear 
theory.  It  is  this  "small-scale  nonlinear  crack  problem"  which  is  studied 
in  detail  in  [  1  ] . 

The  particular  material  considered  in  [l]  belongs  to  a  special  class 
of  incompressible  elastic  materials.  A  characterizing  property  of  this  class 
is  that  each  material  in  it  has  a  strain  energy  density  that  is  completely 
determined  by  the  shear  stress  response  in  simple  shear.  The  specific 
member  of  this  class  used  in  the  pilot  example  of  [  l]  was  chosen  with  two 
requirements  in  mind:  it  was  to  allow  for  a  potential  loss  of  ellipticity, 
and  it  was  to  permit  a  fully  explicit  solution  of  the  small-scale  nonlinear 
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crack  problem  described  above. 

The  anti-plane  shear  field  near  the  tip  of  a  crack  had  been  analyzed 
earlier  in  [2]  for  a  class  of  incompressible  elastic  materials  which  remain 
elliptic  regardless  of  the  severity  of  the  deformation.  For  such  materials, 
the  elastostatic  field  is  infinitely  smooth  in  the  interior  of  the  body,  although 
the  displacement  gradients  as  well  as  at  least  some  of  the  stresses  are 
unbounded  at  the  crack-tip.  In  contrast,  the  example  discussed  in  [l] 
shows  that  a  material  capable  of  losing  ellipticity  will  in  general  give  rise 
to  a  field  near  a  crack-tip  in  which  even  the  first  derivatives  of  displacement 
are  no  longer  continuous  everywhere  in  the  interior.  In  fact,  the  principal 
distinguishing  feature  of  the  results  in  [l]  is  the  appearance  in  the  elasto¬ 
static  field  of  two  curves,  issuing  from  the  crack-tip,  symmetrically 
located  with  respect  to  the  crack-axis,  and  terminating  in  the  interior  of 
the  body,  across  which  stresses  and  displacement  gradients  suffer  jump 
discontinuities.  Some  of  the  local  properties  of  such  "elastostatic  shocks" 
have  been  discussed  elsewhere  [3,4], 

In  the  present  paper,  we  treat  the  small-scale  nonlinear  crack 
problem  for  a  class  of  incompressible  materials  capable  of  losing  ellipticity  — 
rather  than  for  a  single  specific  material.  Although  the  analysis  is  based, 
as  in  [  l],  on  a  hodograph  transformation,  the  emphasis  here  is  on  qualitative 
results  and,  in  contrast  to  the  arguments  used  in  [l],  the  methods  em¬ 
ployed  do  not  depend  on  explicit  solvability  of  the  equations  involved. 

1 .  Special  classes  of  incompressible  elastic  materials 

A  homogeneous,  isotropic,  incompressible  elastic  material  is 
characterized  by  a  strain  energy  density  (or  elastic  potential)  W  which  is 
in  general  a  function  of  the  first  two  of  the  three  fundamental  scalar 

I 
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invariants  I  ,I_,I  of  the  left  (or  right)  Cauchy- Green  deformation  tensor: 

1  b  J 

W  =  Wflj.I^).  Since  only  locally  volume-preserving  deformations  are 
admissible  for  such  materials,  one  always  has  I^  =  1.  In  the  undeformed 
configuration,  I ^  =  1^=  3 ;  in  any  locally  volume -preserving  deformation, 

Ij*3  and  I2^3. 

The  special  class  of  incompressible  materials  with  which  we  shall 
be  exclusively  concerned  consists  of  those  for  which  W  is  independent  of  I^: 

W  =  W(I1)  ,  IjS3  ;  W(3)  =  0  .  (1.1) 

At  infinitesimal  deformations,  the  shear  modulus  for  such  a  material  is 
given  by 

U  =  2W,(3)>0  ,  (1.2) 

where  the  prime  indicates  differentiation.  One  can  show  that  the  Baker- 
Ericksen  inequality  reduces  to 

W,(I1)>0  for  Ij  >  3  (1.3) 

for  materials  of  this  kind. 

Suppose  that  a  body  is  composed  of  a  material  with  the  strain  energy 
density  (1.  1),  and  assume  that  the  region  occupied  by  the  body  in  its  un¬ 
deformed  state  is  a  cube.  Choose  a  rectangular  cartesian  coordinate  frame 
with  its  origin  at  one  vertex  of  the  cube  and  axes  collinear  with  three  of 
the  edges.  Consider  the  simple  shear  in  which  a  particle  with  coordinates 
Xj,X2,Xj  in  the  undeformed  state  is  carried  to  the  point  with  coordinates 

yry2'y3'  where 


: 
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yl  =  xl'  y2=X2'  y3  =  x3+kx2  ; 


(1.4) 


the  nonnegative  constant  k  is  the  amount  of  shear.  *  In  simple  shear^ 


the  invariant  Ij  is  given  by 


Ij=3  +  k  . 


(1.5) 


and  the  shear  stress  r^2=T23  *s  re^ate^  to  ^e  amount  of  shear  by 


T32=T23=T(k,"2kW/(3  +  k  1  ,  k2  0  . 


(1.6) 


The  re sponse  curve  in  simple  shear  is  the  graph  of  T(k).  F rom  ( 1 .  1 ), 
(1.5)  and  (1.6)  it  is  clear  that  W(Ij)  is  determined  by  T(k)  according 
to  the  formula 


/V3 

W(I1)=  j  T(k)dk  ,  1^3  . 


(1.7) 


The  simplest  case  for  which  (1.  1)  holds  is  that  of  the  Neo-Hookean 


material: 


W(I1)=  |(Ir3)  , 


(1.8) 


for  which 


T(k)=uk  ,  k^O 


(1.9) 


The  deformation  (1.4)  is  of  course  locally  volume-preserving. 
‘See  Section  2  of  [2], 


j  ■ 
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A  second  example  is  furnished  by  the  "power-law  materials"  introduced  in 
[  2 ] ;  here 

W(Ii)afe{ri+s(Ir3)’,n-1;  •  !iS3  *  d.10) 

where  u,  n,  b  are  positive  material  constants.  The  response  curve  in 
simple  shear  in  this  case  is  described  by 

x(k)  =u(l +^k2^  k  ,  k^O  .  (1.11) 

Figure  1  illustrates  the  response  in  shear  for  various  values  of  the  "harden¬ 
ing  parameter"  n.  A  power-law  material  hardens  or  softens  in  shear  ac¬ 
cording  as  n>l  or  n<l.  The  case  n=l  corresponds  to  the  Neo-Hookean 
material. 

The  special  material  considered  in  [l]  has  a  strain  energy  density 
given  by 


where  u  is  a  positive  constant.  The  corresponding  response  in  simple 
shear  is: 


Figure  2  displays  the  graph  of  (1.  13).  Note  that  T'(k)  suffers  a  jump 
discontinuity  at  the  peak  corresponding  to  k  =  1  (or  1^  =  4),  and  observe 
as  well  that  T(k)-*0  as  k-»oo. 

In  the  present  paper  we  shall  be  concerned  with  a  subclass  of  the 
materials  characterized  by  (1.  1).  This  subclass  contains  all  materials 
satisfying  the  following  conditions: 

(i)  The  strain  energy  density  W  satisfies  (1.  1),  (1.  2)  and  (1. 3), 
and  the  shear-stress  response  T(k)  given  by  (1.6)  is  three 
times  continuously  differentiable  for  k^O. 

(ii)  T'(k)>0  for  0£k<l  ,  t'(1)  =  0  , 

(1.14)1 

T'(k)<0  for  k>  1  ,  t "( 1 ) < 0  . 

(iii)  -T(k)<  kTr(k)<  t(k)  ,  k>0  .  (1.15) 

(iv)  As  k-*oo  f 

T(k)~Too  +  hk‘1 +e  ,  (1.16) 

where  the  constants  t  ,  h  and  €  are  such  that 

CD 

,  h>0  ,  0<  e  <  1  .  (1.17)2 

According  to  (i),  we  have  T(k)~|jk  as  k-»0.  Condition  (ii) 
asserts  that  the  response  curve  in  simple  shear  rises  monotonically  to 

*  Condition  (i)  and  the  first  three  of  the  four  requirements  in  condition  (ii) 
automatically  imply  t"(1)£0.  We  exclude  the  possibility  t"(1)  =  0  only 
for  reasons  of  simplicity. 

2 

It  may  be  confirmed  that  (1.16),  (1.17)  are  consistent  with  (1.  15). 
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a  maximum  at  k=  1,  then  declines  to  the  limiting  value  .  We  set 

Tm=T(l)  ;  (1.18) 

evidently  is  the  maximum  stress  attained  in  simple  shear. 

The  right  half  of  the  double  inequality  in  condition  (iii)  holds  trivially 
for  1  by  condition  (ii),  but  for  k<  1  it  represents  an  additional  restriction 
and  asserts  that  the  "shear  modulus"  T(k)/k  at  an  amount  of  shear  k  is  a 
decreasing  function  of  k.  This  implies  that 

Tmk^  T(k)Suk  ,  Ofiks  1  ,  (i.  19) 

which  in  turn  requires  that 

.  (1.20) 

The  left  half  of  the  inequality  (1.  15)  follows  from  (1.  14)  for  0<k£  1, 
but  not  for  k>l.  It  requires  that  kT(k)  be  increasing  and  hence  that 

T 

r(k)>-jP  for  k>  1  .  (1.21) 

Finally  condition  (iv)  permits  the  limiting  value  of  T(k)  as  k-oo 
to  be  either  positive  or  zero  and  restricts  the  rate  of  approach  to  this  limit. 

The  power-law  materials  described  by  (1.11)  satisfy  conditions  (i)- 
(iv)  provided  0<n<l/2,  if  one  takes  b=n/l  -  2n,  so  that  the  maximum  in 
the  response  curve  is  at  k=  1.  In  this  case  one  has  Too=0  ,  h  =  u(l  -  2n)*~n, 

€  =  2n  and  Tm  =  b[  2(1  -  n)/(l  -  2n)]n_  1  . 

*  The  location  of  the  maximum  at  k  =  1  is  a  matter  of  convenience  only. 

ft 
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The  material  on  which  the  example  in  [l]  is  based,  although  broadly 
similar  to  materials  in  the  subclass  described  by  conditions  (i)  -  (iv)  in 
that  its  shear  stress  response  curve  rises  to  a  peak  at  k  =  1  and  then  declines, 
does  not  in  fact  belong  to  that  subclass.  It  fails  to  be  smooth  enough  to 
satisfy  condition  (i),  and  violates  the  strict  inequality  on  the  right  in  condition 
(iii)  for  0<  k<  1 . 

2.  The  small-scale  nonlinear  crack  problem 

Let  SI  be  the  three-dimensional  domain  consisting  of  the  exterior 
of  a  plane  infinite  strip  of  width  2c.  Introduce  cartesian  coordinates 
Xj,  x^,  x  in  such  a  way  that  this  strip,  which  represents  a  crack,  lies 
in  the  plane  x^=0,  the  x^-axis  being  parallel  to  the  edges  of  the  crack 
and  the  origin  midway  between  the  edges.  The  cross-section  £  of  9 
in  the  plane  x^=0  is  evidently  the  exterior  of  the  segment  -c^x^  c  of 
the  x  j  -  axi  s . 

Suppose  that  the  interior  of  a  body  of  incompressible  material 
characterized  by  (1.  1)  occupies  the  domain  ft  in  an  undeformed  configuration, 
and  suppose  further  that,  in  the  absence  of  body  forces,  the  solid  is  subjected 
at  infinity  to  a  simple  shear  parallel  to  the  edges  of  the  crack.  1  The  dis¬ 
placement  components  u^  must  thus  satisfy^ 

ua=°0),  u3=kx2+o(l)  as  x  x  -oo  ,  (2.1) 

where  the  nonnegative  constant  k  is  the  amount  of  applied  shear.  The 

This  is  the  loading  case  known  as  "Mode  III"  in  the  terminology  of  fracture 
mechanics . 

2 

Latin  and  Greek  subscripts  have  the  respective  ranges  (1,2,3)  and  (1,2), 
repeated  subscripts  are  summed,  and  a  subscript  preceded  by  a  comma 
indicates  differentiation  with  respect  to  the  corresponding  x-coordinate . 
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deformed  faces  of  the  crack  are  to  remain  free  of  traction. 

For  materials  governed  by  the  strain  energy  density  (1.  1),  it  is 
shown  in  [2]  that  the  field  equations  and  boundary  conditions  associated 
with  the  crack  problem  described  above  are  consistent  with  the  assumption 
that  u^=  0  and  u3=u(xj,x.,)  on  3  ,  corresponding  to  a  state  of  anti -plane 
shear  with  out-of-plane  displacement  ufx^x^).  Moreover,  the  three- 
dimensional  boundary-value  problem  is  reducible,  as  shown  in  [2],  to  the 
following  two-dimensional  problem  for  u  on  the  cross-sectional  domain  £ 


[W'dj^],^0  on  £  , 


(2.2) 


with 


I1=  3  +  |vu|  2  ,  |7u|2=u^u>q  , 


(2.3) 


and 


u(x, ,  x,)  =  kx,+ o(l)  as  x  x  -00  , 

\  c  £.  a  a 


(2.4) 


u  2(xi»0±)  =  0  for  -c<Xj<c  . 


(2.  5) 


It  is  also  shown  in  [2]  that  the  components  of  true  stress  t_  satisfy 


T3<rTa3=2W'(Il)u,„  •  . 


(2.6) 


ToP=0  • 


(2.7) 


with  Ij  given  by  (2.3). 

When  the  Bake r- Ericksen  inequality  holds,  as  we  assume  here. 


i/  W  -  v  v 

/ 
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it  is  possible  to  show  that  the  differential  equation  (2.  2)  is  elliptic  at  a 
solution  u  and  at  a  point  (xyX^)  if  and  only  if 

W'(3  +  k2)  +  2k2W"(3  +k2)>0  ,  k  =  |7u(Xj,  x^  |  .  (2.8) 

In  view  of  (1.6),  this  condition  for  ellipticity  is  equivalent  to  the  assertion 
that  the  response  curve  in  simple  shear  must  have  positive  slope  at  an 
amount  of  shear  equal  to  the -magnitude  of  the  local  displacement  gradient 
^ufxj.x^).  It  follows  that  (2.2)  can  never  suffer  a  loss  of  ellipticity  if  the 
response  function  T(k)  for  simple  shear  is  monotone  strictly  increasing. 

On  the  other  hand,  for  the  class  of  materials  satisfying  conditions  (i)  -  (iv) 
at  the  end  of  the  preceding  section,  (2.  2)  will  lose  ellipticity  whenever 
!^u(Xj,x7)!  SI.  In  particular,  for  the  power-law  materials  characterized 
by  (1.  10),  (1.  11),  a  loss  of  ellipticity  is  possible  only  if  n<  1/2. 

In  the  presence  of  a  potential  loss  of  ellipticity,  it  is  natural  to 
seek  the  solution  of  the  boundary-value  problem  (2.2)  -  (2.5)  in  a  class 
of  functions  whose  smoothness  is  less  than  that  which  one  would  otherwise 
naturally  demand.  In  fact  we  shall  merely  require  that  u  be  continuous 
and  have  piecewise  continuous  first  and  second  partial  derivatives  on  £. 

In  addition  u  is  to  be  bounded  inside  any  circle  centered  at  a  crack-tip. 
Finally,  the  limits  ?u(xj,0±)  must  exist  and  be  continuous  for  -c<Xj<c. 

The  preceding  smoothness  requirements  allow  for  the  possibility 
of  finite  jump  discontinuities  in  Vu—  and  hence  also  in  the  stresses  — 
across  curves  in  £.  It  is  necessary  to  stipulate  that  across  such  a  curve 
J  ,  W'(Ij)8u/8n  shall  be  continuous*,  where  1^  is  given  by  (2.3)  and 

*This  is  equivalent  to  the  requirement  that  the  tractions  be  continuous  across 
J  ,  which  in  turn  is  demanded  by  equilibrium.  See  Section  1  of  [  1  ]  and 
Section  3  of  [3  ] . 
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9u/9n  is  a  derivative  of  u  normal  to  J . 

For  small  values  of  the  amount  of  applied  shear  k  ,  one  would 
expect  that,  away  from  the  crack-tips,  the  solution  of  the  foregoing  boundary- 
value  problem  would  be  well- approximated  by  the  solution  of  the  corresponding 
linearized  problem.  1  In  fact,  arguments  sketched  in  [l]  and  [2]2  suggest 
that,  for  small  k,  an  approximation  to  the  solution  of  (2.2)  -  (2.5)  which  is 
uniformly  valid  near  the  right  crack-tip  is  obtained  by  replacing  the  domain 
£  with  the  domain  7)  exterior  to  the  nonpositive  Xj-axis  (see  Fig.  3),  re- 
quiring  (2.2)  to  hold  on  £  ,  enforcing  the  free  surface  condition  (2.5)  on 
the  semi-infinite  crack  -oo<x^<0  ,  x^=  0  ,  and  matching  u  and  its  first 
derivatives  at  infinity  to  the  corresponding  near-field  quantities  calculated 
from  the  exact  solution  of  the  original  problem  according  to  the  linearized 
theory.  This  leads  to  the  small-scale  nonlinear  crack  problem  in  which 
one  seeks  a  solution  u  of  (2.2)  on  7,  subject  to  the  boundary  conditions 

u  ,  Xj<0  (2.9) 


and  the  matching  conditions 


.  „  .1/2  .  9 

u~k(2cr)  sin 


A 


u  j~-kc(2cr)  *^sin  ,  u  2~kc(2cr)  ^^cos-|  ,  ^  (2.10) 


as  r  -*oo  ,  -ir£9s-rr  . 


Upon  linearization,  the  boundary-value  problem  (2.2)-(2.5)  passes  over  to 
the  elementary  problem  for  Laplace's  equation  corresponding  to  irrotational 
flow  of  an  inviscid  incompressible  fluid  past  a  flat  plate  at  a  90°  angle  of 
attack. 

See  also  Rice  [  5] . 
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Here  r,  9  are  polar  coordinates  at  the  crack-tip  (Fig.  3),  so  that 

Xj=rcos9  ,  x^=  r  sin  6  ,  r^O  ,  -its  8^tt  .  (2.11) 

Furthermore,  u  is  to  satisfy  relaxed  regularity  conditions  strictly  analogous 
to  those  imposed  in  connection  with  the  original  crack  problem. 

When  referred  to  the  dimensionless  variables 

xa=xa/ck^  ,  r  =  r/ck^  ,  u  =  u/ck^  ,  (2.12) 

the  boundary -value  problem  (2.2),  (2.9),  (2.  10)  assumes  the  following  non- 
dimensional  form  in  which  k  does  not  appear  explicitly: 

[  W'(3  +  |Vu  j  2)u  ^=0  on  £  ,  (2.13)1 

u  t  Xj<0  ,  (2.14) 

_  ,,-,1/2  .  6 

u~(2r)  sin  , 

/->  v-1  /2  .  9  ,,  ,-1/2  e 

u  j - (2r)  sm  2  »  u  2~'2r  cos  2  • 

as  r  -cd 

The  foregoing  small-scale  nonlinear  crack  problem  was  solved  in 
[2]  for  the  limiting  elliptic  case  n=  1/2  of  a  power-law  material  governed 
by  (1.  10),  (1.11),  (see  Fig.  1),  and  in  [  1  ]  for  the  :  laterial  described  by 
(1.12),  (1.13),  (Fig. 2). 

*  Subscripts  preceded  by  a  comma  are  now  understood  to  indicate  partial 
derivatives  with  respect  to  the  dimensionless  coordinates  x  ,  and  V  stands 
for  the  gradient  operator  with  components  9/9xa  • 
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3 .  An  algorithm  for  the  construction  of  solutions  to  the  displacement 
equation  of  equilibrium. 

A  formal  algorithm  for  the  construction  of  a  solution  to  the  small- 
scale  nonlinear  crack  problem  (2.  1 3 ) - ( 2 .  15)  can  be  generated  with  the  help 
of  a  hodograph  transformation,  in  which  the  first  derivatives  of  u  are 
introduced  as  new  independent  variables,  and  one  is  led  to  an  explicitly 
solvable  linear  boundary-value  problem  in  the  hodograph  plane  for  the 
Legendre  transform  of  u  .  For  materials  which  remain  elliptic  at  all 
deformations,  the  mapping  between  the  physical  and  hodograph  planes 
is  one-to-one,  and  the  formal  procedure  indeed  produces  an  exact  solution 
which  is  infinitely  smooth  on  £>  .  This  analysis  is  carried  out  in  detail 
in  [2]  for  the  limiting  elliptic  case  n  =  1/2  of  the  power-law  material 
(1.  10). 

In  [  1  ]  the  formal  hodograph  scheme  is  applied  to  the  boundary- value 
problem  (2.  13)-(2.  15)  for  the  case  of  a  material  defined  through  (1.  12). 

In  this  instance,  a  loss  of  ellipticity  occurs  and  is  accompanied  by  a  break¬ 
down  in  the  inve rtibility  of  the  hodograph  transformation.  As  shown  in  [  l], 
even  in  these  circumstances  the  formal  hodograph  procedure  can  be  used 
to  find  solutions  of  the  differential  equation  (2.  13).  These  can  in  turn  be 
combined  to  construct  a  solution  of  relaxed  smoothness  to  the  boundary- 
value  problem. 

We  shall  not  repeat  the  description  of  the  hodograph  transformation 
here;  instead  we  merely  cite  from  [l]  the  resulting  formal  algorithm  and 
analyze  it  for  incompressible  materials  governed  by  (1.  1)  and  satisfying 
conditions  (i)  -  (iv)  stated  near  the  end  of  Section  1.  Let  T(k)  be  the 
shear-stress  response  function  for  such  a  material,  and  define  functions 

i 
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P  and  Q  by 


P(p)=  f  ,  Q(p)=-pP'(o)  ,  p>0  ,  (3.1)1 

;  t2?(t) 


where 


T(t)  =  r(t)/u  .  (3.2) 

According  to  the  hodograph  scheme,  the  formal  solution  u  of  the  boundary- 
value  problem  (2.  1 3 ) -  (2 .  15)  is  given  implicitly  by 

u(xrx2)  =  -pQ(p)cosco  ,  r>0  ,  ,  (3.3) 


where  for  given  (Xj.x^)  or  (r,  9)  ,  p  and  cp  are  to  be  determined  from 


x j  =  r  cos  6  =  P(p )  -  Q(P )cos  cp  , 

X2=irsin9  =  -Q(p)sincpcos  o  ,  p>0  ,  . 


>  (3.4) 


It  can  be  shown  that  for  an  elliptic  material,  (3.4)  are  uniquely 
invertible:  for  each  r>0  and  each  6  in  [  -tt , it ] ,  there  are  unique  values 
p=p(r,  9),  cp  =  cp(r,  9)  with  p>0  and  O^cpSir  such  that  (3 . 4)  hold.  Substitu¬ 
tion  of  these  values  of  p ,  co  into  (3.3)  furnishes  the  solution  u  . 

For  the  materials  of  interest  here,  (3.4)  will  in  general  fail  to  be 
uniquely  invertible  for  certain  values  of  (xj,x2)-  The  first  step  in  con¬ 
structing  a  solution  to  the  small-scale  nonlinear  crack  problem  (2.  13)- (2.  15) 


*Note  that  (1.  16)  assures  the  convergence  of  the  integral  in  the  definition  P. 
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consists  in  determining  the  number  of  roots  (p,;p)  of  (3.4)  for  each  value 
of  (Xj.xJ. 

A  geometrical  version  of  this  issue  is  arrived  at  by  observing  that 
(3.4)  may  be  written  in  the  equivalent  alternative  form 


whe  re 


Xj=  r  cos  9  =  a(p)  -  b(p  tcos  2cp 
X2  =  'rsin9  =-b(p)sin2cp  , 


y  (3.5) 


OscpSir  ,  p>0  ,  I 


a(p)  =  P(p)  -  -|q(d)  ,  b(p)  =  -|Q(p)  .  (3.6) 

F rom  (3.1),  (3.2),  (1.6)  and  (1.3)  it  follows  that  b(p ) >  0  for  p  > 0  ,  so 
that  (3.5)  represents,  for  each  fixed  p>0  ,  the  parametric  equations  of 
a  circle  T(o)  centered  at  Xj=a(c)  ,  x^  =  0  with  radius  b(p).  The  cartesian 
equation  for  this  circle  is 

r(p):  [xr  a(p)]2+X2  =  b2(p)  .  (3.7) 

Thus  for  given  values  of  Xj  and  x^  ,  (3.4)  has  a  solution  (p.,.,  CD.J  with 
P*>0  »  O^cp^Sir  ,  if  and  only  if  the  circle  HP*)  passes  through  the 
point  (xj.x^).  Moreover  if  (p^cp^)  and  (p^.,  cp^.,.)  satisfy  (3 . 4)  for  given 
(xrx2)  ,  then  cp.,.  =  cp;,;j,.  .  Thus  the  number  of  admissible  roots  ( P ^ » CP*. ) 
of  (3.4)  coincides  with  the  number  of  circles  T(p)  which  pass  through  the 
given  point  (x^.x^)  or,  equivalently,  with  the  number  of  roots  p^  of 
(3.7). 


-17- 


Let  be  a  subdomain  of  £i  on  which  is  a  twice  continuously' 
differentiable  positive  function  of  (Xj.x^)  satisfying  (3.7).  Note  that  for 
such  a  p.;.  , 

a(pj  -  b(pJ.)S  x a(p„.)  +  b(p„.)  ,  (3.8) 

since  the  left  and  right  members  of  these  inequalities  are  the  abscissas 
of  the  intercepts  of  T’fp.J  with  the  x^-axis.  Define  cd...  through 

cos  cpj.=  -sgn  x^ 

One  shows  easily  that  o;|c,cp.,.  satisfy  (3.4)  on  .  Further,  if  c p.,,  is 
twice  continuously  differentiable  on  ,  the  function  U  defined  on  by 

u(xj,x2)=  -p  ...(x^ ,  x^JQfo  ...(x  ^  ,  j^Dcos  cp.:.(Xj ,  x^)  (3.10) 

can  be  shown  by  direct  calculation  to  be  a  solution  of  the  partial  differential 
equation  (2.  13)  on  . 

In  the  next  section  we  investigate  the  roots  p...  of  (3.7)  and  their 
domains  of  definition  . 

4.  Geometry'  of  the  circles  T(p) 

In  order  to  investigate  the  roots  p  of  (3.7),  we  need  some  properties 
of  the  functions  a(p)  and  b(o)  which,  according  to  (3.7),  determine  the 
circles  I"(p). 

From  the  definitions  (3.6)  and  (3.  1)  it  follows  that 
^sgnx^  is  one  if  *  zero  if  an<*  minus  one  ^ 


a(p.J  +  b(p . 


-  11/t 


-  X. 


2b(pJ 


0  S  CDj.  ^  TT 


(3.9) 
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oo 

=  \  ~r^ - '  b(p)  =  T^”T  •  p>0  • 

D  t  T(t)  2pt(g)  2dt(p) 


a(p) -  ]  2_ 


1 


1 


(4.  1) 


,(p)=p iM-M  ,  b'(D)=_£iM+iM  ,  0>0  , 

2p  t  (p) 


,  2-2/  > 
2p  t  (p) 


(4.2) 


and  hence  in  particular  that 


a'(p)  -  b'(p)  =  — ,  P  >  0 
PT  (P) 


(4.3) 


Invoking  (1.15),  one  finds  from  (4.2)  that 


a'(p ) <  0  ,  b'(p)<0 

,  p  >0  , 

(4.4) 

while  (4.3),  (1.  14)  lead  to 

r 

>0  , 

f-H 

V 

Cl 

V 

o 

a'(p)  -  b'(P  )•< 

=  0  , 

p  =  i  . 

(4.  5) 

<0  , 

p  >  i  . 

We  shall  also  require  the  asumptotic  properties  of  a(p)  ,  b(p)  for 

large  and  small  p.  The  asymptotics  for  large  k  of  t ( k)  given  by  (1.  16), 

(1.  17)  lead  via  (4.  1)  to  the  following  estimates  for  a(p),  b(p): 

If  t  >0:  as  p  -ao  , 
ao 


a  (d  ) 


U  2  _  uhe _ 1 

2t  p  ,2  .  2  -  e  ' 

oo  v  2t00(2  -  e)  p 


b(p) 


u  1  uh  1 


2t  p  ,2  2-e 

oo  2t  p 
oo  y 


>  (4.6) 


J 
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Lf  t  =0:  as  p  -  oo 

CD 


,  ,  u(2  -  e)  1  u  1 

a(?)  Zhl  c  '  b  P  ~  2h  ~ 

P  P 


(4.7) 


From  (1.6),  (1.2),  (3.2)  one  finds  the  behavior  of  t(k)  for  small 


k  to  be  given  by 


F(k)~  k-dk  as  k-0  , 


d  =  -2W"(3)/u  . 


(4.8) 


Ruling  out  a  degeneracy,  we  shall  assume 


Using  (4.  8),  (4.  9)  in  (4.  1)  leads  to 


(4.  9) 


i(o)  ~  -d  log  2  ,  b(p ) - ~  as  D-0  + 

2o^ 


(4.  10) 


From  (4.6),  (4.7)  it  may  be  noted  that  in  either  of  the  cases  t  >0 


or  r  =  0  ,  one  has 
CD  ' 


a(p)-0  ,  b(p)-*0  as  c-co  , 


(4.  11) 


while  from  (4.  10)  it  follows  that 


a(p)-»  +  oo  ,  b(p)-  +  oo  as  P-0  + 


(4.  12) 


On  the  other  hand,  (4.  10)  also  shows  that 


The  fact  that  d  ^  0  is  a  consequence  of  (1 .  15). 


/  k *■;  l  - 
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a(p )  -  b(p )  -  -  oo  as  p-*0+  .  (4.13) 

From  (4.4),  (4.  12)  and  (4.  11)  one  infers  that  a(p)  and  b(p)  decrease 
monotonically  from  +oo  at  p  =  0+  to  zero  at  p  =+oo.  In  particular,  both 
functions  are  positive  for  0<p<oo. 

According  to  (4.  13),  (4.  5)  and  (4.  11),  a(p)  -  b(c)  increases  mono¬ 
tonically  from  -  oo  at  p  =  0+  to  a  maximum  at  p  =  1  and  then  decreases 
steadily  to  zero  at  p=+oo  (Fig.  4).  Since  the  maximum  value  of  a(o)  -  b(c) 
is  necessarily  positive,  one  has 

a( 1 )  -  b( 1 ) > 0  .  (4.  14) 

Finally  it  is  clear  that  there  is  a  unique  number  Pq  such  that  0  <  P  0  <  i  and 

a(p0)  -  b(cQ)  =0  .  (4.15) 

The  circles  T(p)  of  (3.7)  are  centered  at  Xj  =  a(p)  ,  x^- 0  >  have 

radius  b(p)  and  intersect  the  Xj-axis  at  Xj=a(p)-b(c)  and  at 

x j  =  a(p )  +  b(p ) .  As  P  increases  from  zero  to  one,  the  center  and  right 

Xj-intercept  of  T(p)  more  from  Xj=  +  oo  to  the  left,  while  the  left 

Xj-intercept  moves  to  the  right.  By  (4.  15),  the  left  intercept  of  "(Pq) 

is  at  the  origin.  The  circle  T(l),  which  plays  a  major  role  in  the  analysis 

to  follow,  is  centered  at  x,  =  a(l)  ,  x,  =  0  ,  has  radius  b(l)=p/2T  ,  and 

i  i  m 

intersects  the  Xj-axis  at  positive  values  of  x^. 

As  p  continues  to  increase  from  unity,  the  center  and  right 
Xj-intercept  continue  to  move  to  the  left,  and  the  radius  continues  to 
decrease,  but  the  left  intercept  Xj=a(c)-b(c)  now  moves  to  the  left  as 
well  (Fig.  5).  For  very  large  values  of  p  ,  T(p)  is  centered  very  slightly 

< 

I 
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to  the  right  of  the  origin  and  has  a  small  radius. 

Two  circles  T(p)  corresponding  to  values  of  p  both  of  which  are 
less  than  one  do  not  intersect,  but  this  is  not  necessarily  the  case  when 
both  values  of  p  are  greater  than  one.  It  is  thus  not  surprising  to  find 
that  the  subfamily  of  circles  T(p)  for  which  p^  1  possesses  an  envelope. 
This  envelope,  which  is  of  decisive  importance  in  the  study  of  the  roots 
p  of  (3.7),  is  a  curve  which  is  tangent  at  each  of  its  points  to  one  of  the 
circles  F(p)  ,  1.  To  find  this  curve  analytically,  we  first  define 

F(p;x1,x2)  =  [xj- a(p)]2+x“  -  b^p)  ,  p>0,  -co<x^<oo  .  (4.16) 

The  envelope  then  consists  of  those  points  (Xj.x^)  for  which  there  is  a 
p  > 0  such  that 

F(p;xJ(x2)  =  0  and  — (p;Xj,x2)=0  .  (4.17) 

Substitution  from  (4.  16)  into  (4.  17)  leads  directly  to  parametric  equations 
for  the  envelope  a £  : 

1/2 

=£:  Xj=ej(p)  ,  x2=±[e2(p)]  ,  p^l  ,  (4.18) 

where  the  functions  e j(i)  ,  e^(p)  are  defined  for  all  p>0  by* 


^Although  only  the  values  of  e^(p)  ,  e^fp)  for  1  enter  into  the  parametric 

equations  for  a£  ,  their  values  for  0<p<l  will  be  needed  later. 
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e  j(o )  = 


a(p)a'(o)  -  b(o  )b'(c ) 
a'(p) 


>  (4.  19) 


e,(p)  =  -^-!£i-,[a'(p)-b'(o)][a'(p)+b'(p))  ,  0<  p  <oo 

2  [  a'(c)] 2 


Only  those  cirdes  T(p)  for  which  0  &  1  participate  in  the  formation  of  the 
envelope  because,  by  (4.4),  (4.5),  e^pJ^O  for  PS1  but  tor 

0<P<1.  Thus  x2  as  given  in  (4.  18)  is  real  only  for  p^l. 

The  envelope  sf  is  symmetric  about  the  x^-axis,  consisting  of  an 
upper  and  a  lower  branch  as  indicated  by  the  sign  alternative  in  (4.  18). 
From  (4.4),  (4.5),  (4.  19)  and  (4.  14)  one  finds  that  e  j  (o )  >  0  for  p^l, 
so  that  =L  lies  in  the  half-plane^  Xj  ^  0.  Moreover,  e  j(  1)  =  a(  1)  -  b(  1)  and 
e 2(  1 )  =  0  ,  so  that  each  branch  of  af  "begins"  when  p  =  l  at  the  left 
Xj-intercept  of  T(l).  A  more  detailed  inspection  of  e  ^ ( p )  ,  near 

p  =  1  shows  that  s t  has  a  vertical  tangent  at  p  =  1  and  its  radius  of  curva¬ 
ture  there  is  b(l).  Thus  *£  osculates  with  F(l)  at  the  left  x^-intercept 
of  this  circle . 


As  P  -*  +  oo  ,  the  two  points  on  sf  which  correspond  to  p  tend  to 
the  origin,  as  one  sees  easily  from  the  asymptotic  behavior  of  e^(p)  and 
e^fc)  for  large  p  deducible  from  (4.6)  or  (4.7),  according  as  T00>®  or 
7^=0.  In  fact,  it  can  be  shown  that  the  upper  and  lower  branches  of  si 
approach  the  origin  with  vertical  tangents  as  p-oo  ,  provided  t  >0. 

If  T(X)=  0  ,  the  upper  and  lower  branches  of  =£  are  asymptotically  tangent 
as  p-  +  oo  to  the  lines  ±  e  ( 1  -  e )"  *  ^ “Sc j  ,  where,  according  to 
(1.  16),  c  is  the  material  parameter  controlling  the  rate  of  descent  of  the 
shear-stress  response  curve  at  large  amounts  of  shear. 


1 


The  origin  (which  would 


"correspond"  to 


p  =oo)  is  taken  to  be  a  point  of  st  . 
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From  (4.  19)  one  finds  that 


e'^p)  =  a'(p)c(p)  ,  e'2(p)  =  2b(o)b'(p)c(p)  ,  p>0  , 

where 

,  fb'\Z,  bb'a"  bb" 
c  =  1  -  — /  + - 5 - 2  • 

a  (a')  (a') 


(4. 20) 


(4. 21) 


By  (4.18),  (4.20),  the  slope  dx^/dxj  of  sf  therefore  satisfies 


x 


dx„ 


2  dx, 


b(p)b'(c)  >n 

a'(p) 


for  eve  ry 


P 


such  that  c(p  )  /  0  . 


(4.22) 


Thus  on  the  upper  (lower)  branch  of  dL  ,  5c ^  is  a  monotone  strictly  increasing 
(decreasing)  function  of  5c^  on  any  arc  of  =£  which  includes  only  points  for 
which  c(p)^0.  At  values  of  o  for  which  c(p)  =  0,  st  has  a  cusp. 

By  (4.21),  (4.  1),  (4.  5),  (4.3),  (1.18)  and  (3.2),  one  has 


c(l) 


[b'uij2 


[a"(X)-b"(ll]  . 

m 


so  that  by  the  last  of  (1.  14), 


(4.22) 


c(l  )<  0  . 


(4.  23) 


It  then  follows  that  c(p)<0  when  p  is  near  unity,  and  hence,  by  (4.  18), 
(4.20),  (4.4),  5cj  at  first  increases  along  dL  as  p  increases  from  p  =  l. 
But  as  o-*oo  ,  the  corresponding  points  on  s£  tend  to  the  origin,  and  thus 
Xj  is  ultimately  decreasing.  One  concludes  that  c(p)  must  change  sign 
at  least  once  for  l<p<oo.  If  the  material  is  such  that  for  all  p  >0  , 
only  one  such  change  of  sign  occurs,  then  there  is  a  unique  number  PC>1 
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such  that 

c(p )  <  0  ,  0<p<Cc  ;  c(dc)  =  0  ;  c(p)>0  ,  P>Pc  .  (4.24)1 

When  (4.24)  holds,  one  may  infer  from  (4.20)  that  ej(p)  and  e.,(p) 

increase  monotonically  with  0  for  0<0<p  ,  attain  maximum  values  at 

c 

p=p  ,  and  thereafter  decrease  monotonically  with  p.  The  asymptotic 
c 

behavior  of  e^(p)  ,  e^p)  ,  which  can  be  determined  from  (4.6)  or  (4.7), 

(4.  10)  and  (4.  19),  is  such  that  e^c)  and  e?(p)  both  tend  to  -  oo  as 
0  -0+  ,  and  both  tend  to  zero  as  p-»  +  oo.  The  zero  of  e^(o)  lies  between 
p  =  Cq  and  p  =  1.  The  graphs  of  e  ^(p)  and  e-,(p)  are  thus  as  sketched  in 
Fig. 6. 

In  the  analysis  to  follow,  it  will  be  assumed  that  the  material  under 
consideration  satisfies  the  "one-cusp  condition"  (4.24).  At  the  end  of 
this  section,  we  shall  discuss  the  interpretation  of  this  condition  in  terms 
of  the  shear-stress  response  curve  of  the  material.  In  the  final  section, 
we  shall  comment  briefly  on  the  effect  of  relaxing  the  restriction  (4.24). 

Suppose  (4.24)  holds,  and  let  DfpJ^O  be  the  distance  from  the 
center  of  the  circle  T(l)  to  the  points  on  JL  corresponding  to  p.  A 
direct  calculation  gives 

D(p  )D'(P  )  =  a'(p )[  a(p )  -  a(  1 )]  c(P )  ,  p^l  ,  (4.25) 

so  that  D'(p)<0  for  l<p<p  ,  D'(p  )  =  0  ,  and  D/(d)>0  for  p>p  .  Thus 

c  c  c 

as  p  increases  from  p  =  1  ,  the  corresponding  point  on,  say,  the  upper 

*  The  asymptotic  behavior  of  c(p)  as  p-*0+  guarantees  that  c(p)<0  for 
p  near  zero. 

2 

Hence,  each  branch  of  sf  has  only  one  cusp. 
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branch  of  af  moves  closer  to  the  center  of  f(l)  until  it  reaches  the  cusp, 
where  D(p)  has  a  minimum.  ^  The  point  then  moves  away  from  the  cusp  — 
which  lies  inside  T(l)  -  in  such  a  way  as  to  increase  its  distance  from  the 
center  of  f(l),  eventually  crossing  T(l)  and  approaching  the  origin. 

By  a  similar  argument,  one  can  show  that  all  points  of  af  except 
the  origin  lie  in  the  interior  of  the  circle  r(pg)  through  the  origin. 

When  (4.24)  holds,  it  can  be  shown  that  there  are  no  points  of 
inflection  on  those  portions  of  the  upper  or  lower  branches  of  a£  for  which 
1  £  c  <  Dc  ,  but  the  remaining  part  of  each  branch  has  at  least  one  point  of 
inflection.  It  is  also  possible  to  prove  that  there  are  no  double  points  on 
= t  when  the  one-cusp  condition  holds. 

A  qualitative  sketch  of  sL  is  shown  in  Fig.  7  for  t  >0  (Case  A) 
as  well  as  for  t  =0  (Case  B). 

CD 

One  can  show  that  the  envelope  s£  is  precisely  the  set  of  points 

(Xj.x,)  at  which  the  Jacobian  determinant  of  the  physical  coordinates  with 

respect  to  the  hodograph  coordinates  vanishes.  It  is  therefore  analogous 

to  a  "limit  line"  in  gas  dynamics.2 

We  have  not  found  a  satisfactory  interpretation  of  the  one -cusp 

condition  (4.24).  One  can  show  that  c(p)  vanishes  at  p=o  >1  if  and 

c 

only  if  the  curve  representing  the  square  of  the  shear  stress  t  in  simple 
shear  as  a  function  of  the  square  of  the  amount  of  shear  p  has  a  point  of 
inflection  at  p  =p^  .  The  presence  of  more  than  one  cusp  on  each  branch 
of  the  envelope  £  would  therefore  correspond  to  "terracing"  in  the 

*One  can  show  that  JL  is  tangent  to  T(pc>  at  the  cusps. 

2See  §105  of  [6]. 
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declining  portion  of  the  latter  carve  (c>l). 

It  can  be  verified  that  the  power-law  materials  (1.  10),  (1.  11)  for 
0<n<l/2  satisfy  the  one-cusp  condition  (4.24). 

5 .  Elliptic  and  hyperbolic  solutions  of  the  differential  equation 

We  turn  now  to  the  question  of  the  number  of  roots  o>0  of  (3.7), 
or  equivalently,  of  the  number  of  positive  zeros  of  F  in  (4.  16)  as  a 
function  of  p  . 

First,  suppose  that  x^  =  x^=0  ,  and  observe  from  (4.  16)  that 

F(p;0,  0)  =  [a(o)  +  b(p)][a(p)  -  b(p ) ]  ,  p>0  .  (5.1) 

Since  a(p)+b(p)>0  for  p>0  and  a(p)-b(p)  vanishes  only  at  p^  (see 

(4.  15)),  it  follows  that  the  only  zero  of  F(p;0,0)  is  Dq. 

_2  _2 

Next,  assume  that  x^+x^^O.  From  (4.  11),  (4.16),  one  sees  that 
F(+oo;Xj ,  x^)  =  Xj  +x^  >0  ,  (5.2) 

while  (4.  10),  (4.  16)  yield 

F(0+;Xj,  x^)  =  -oo  .  (5.3) 

Thus  for  any  given  (xj.x^)  ,  there  is  at  least  one  positive  zero  p  of 
FfpjXj.x^).  Note  that  this  result  does  not  depend  on  the  one-cusp  condition 
(4.24). 

In  order  to  investigate  the  number  of  zeros  of  F  ,  it  is  convenient  to 
record  first  some  properties  of  this  function.  From  (4.  16),  (4.  19)  one 
obtains 


/  '5  -  V. 


' — w 
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(cjXji^)  =  -2a'(p)[x1-e1(c)]  ,  P>0,  -oo<xo,<ao  (5.4) 

Thus  by  (4.4),  the  graph  of  F  vs.  P  for  given  (Xj.x^  has  a  horizontal 
tangent  when  p  is  such  that  Xj  =  ej(p).  One  sees  from  Fig.  6, which  depends 
on  (4.24),  that  when  x^  £  0  ,  there  is  precisely  one  such  horizontal  tangent  — 

located  at  a  value  of  p  between  zero  and  one.  If  0<Xj<e^(p  )  ,  two 

horizontal  tangents  occur  at,  say,  P  =Pj  •  &  ~D2  Pj<~Pc<^p2  '  ^r^en 

p^>l  ,  both  Pj  and  p^  correspond  to  points  on  £  which  are  intersections 
of  £  with  the  vertical  line  through  the  point  (Xj,x->).  When  Pj<1  ,  only 
p^  refers  to  a  point  on  £  with  this  property.  *  If  xi  =  ei^cc^  >  is  a 

horizontal  tangent  in  the  graph  of  F  only  at  p  =  c  ,  and  if  x^>Pc  ,  there 

are  no  horizontal  tangents. 

Equations  (5.4),  (4.  16)  and  (4.  19)  show  that,  at  a  horizontal  tangent 

of  F, 


F(p;x1,  x2)  =  x2  -  e2(p) 


dp  ($''xi’x2^  "  0 


F rom  the  graph  of  e2(p)inFig.6,  it  is  clear  that  ^  ®  if  P ^  1  »  so 

that  F(p  ;Xj ,  x2)  >  0  at  a  horizontal  tangent  for  which  C<1.  If  P  s  1  in 
(5.5),  then  F(p;Xj,x2)  is  positive,  zero  or  negative  according  as  the  point 
(Xj,  |x2|)  lies  above,  at  or  below  that  point  on  the  upper  branch  of  the 
envelope  £  which  corresponds  to  p.  The  latter  point  is  a  point  of  inter¬ 
section  of  £  with  the  vertical  line  through  (xj,x2). 

From  (4.  16)  we  observe  that  F(l;xj,x2)  is  positive,  zero  or  negative 
according  as  (Xj,x2)  lies  outside,  on  or  inside  the  circle  T(l).  Note  that 

^Recall  that  £  contains  no  points  corresponding  to  values  of  p  <  1  (see  (4.  18)). 
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fore  ach  fixed  p  and  Xj,  F(p;x1,x2)  is  a  monotone  increasing  function 
of  x2  for  x2>0. 

The  preceding  considerations  lead  to  the  following  conclusions 
concerning  the  dependence  of  the  number  and  location  of  the  zeros  of  F 
on  the  point  (xj.x.,): 


Table  1 

Location  of  (x^,x2) 

Number  of 
ze ros  of  F  in  (0,1) 

Number  of 

zeros  of  F  in  ( 1 ,  oo) 

outside  JL  ,  outside  T(l) 

1 

0 

outside  =d  ,  inside  f(l) 

0 

1 

inside  :£  ,  outside  r(l) 

1 

tm 

inside  d  ,  inside  F(l) 

0 

3 

One  observes  in  particular 

that,  at  any  point  (Xj,x2) 

which  lies  outside  sd 

there  is  a  unique  zero  of  F 

,  while  at  any  point  (Xj, 

x0)  in  the  interior  of 

si  ,  there  are  three  distinct  zeros. 

If,  at  a  given  value  (x^.x?,)  of  (Xj,x2)  ,  p  is  a  zero  of  F  ,  the 
implicit  function  theorem  assures  the  existence  of  a  function  p.,.(xj,x2), 
defined  and  smooth  in  a  neighborhood  of  (x°j,x^),  such  that 
F(p...(xj  (  x2);Xj  ,x2)  =  0  in  this  neighborhood,  unless  3F/9p  vanishes  at 
p  ,  x°j,3c^  .  Thus  differentiability  of  p ^  can  breakdown  at  (xj,x2>  only 
if  (3Ej,3c2)  coincides  with  the  point  on  the  envelope  sL  that  corresponds  to 
0  .  Since  there  are  no  points  on  sd  corresponding  to  values  of  p<l  ,  no  such 
breakdown  in  smoothness  can  occur  if  p<  1. 

We  now  define  two  functions  ,  p2  on  the  (xj ,  x2)-plane  as  follows: 

p  j(Xj ,  x2)  =  the  smallest  zero  of  F(p;x^,x2)  , 


(5.6) 
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P2(Xj,  x^)  =  the  largest  zero  of  Ffpjx^.x^)  . 


(5.7) 


Since  F  has  only  one  zero  for  points  (x^x^  which  lie  strictly  outside 
sc  ,  one  has  p  P£  in  the  exterior  of  s£  .  From  (5.  1),  (4.  15)  one  has 
p  j(0, 0)  =  p ^(0 ,  0)  =  Pq  ;  p  j  is  continuous  at  the  origin,  but  p^  is  not.  In 
fact,  one  has  P^Cxj.OI-^Pq  as  x^-»0-  ,  while  p ^(Xj ,  0)  -  +  co  as  Xj-0+  . 
Also,  ^\=d2  =  ^c  ^  Oxj.x^)  coincides  with  either  cusp,  and  p^  ,  p^ 
are  continuous  at  the  cusps.  One  can  show  that  Pj  suffers  jump  discon¬ 
tinuities  at  all  points  on  that  arc  of  the  envelope  sf  which  joins  the  cusps 
by  way  of  the  left  x^-intercept  of  T(l)  ,  except  for  the  cusps  themselves. 
Except  for  points  on  this  arc  of  £,  pj  is  three  times  continuously  differ¬ 
entiable.  On  the  other  hand,  P^  has  jump  discontinuities  at  all  points 
of  that  branch  of  £  which  joins  the  cusps  via  the  origin,  the  cusps  ex¬ 
cluded,  and  p.,  has  continuous  third  derivatives  elsewhere.  Both  p.  p 
are  even  functions  of  x^  for  each  x^  . 

One  infers  from  the  foregoing  discussion  that 


Pj  =  p2<i  outside  both  af  and  T(l)  , 
p  1  =  D ^ ^  1  outside  but  inside  F(l)  , 


>  (5.8) 


y 


c  i  >  1  ,  P2>1  inside  both  £  and  T(l)  , 
Pq<Pi<1  f  p2>^  inside  sL  but  outside  T(l)  , 


> 


J 


(5.9) 


and 


Pq<Pj<Pc  ,  P2>PC  inside  ^  . 


(5.  10) 
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As  for  fixed  Xj  ,  one  shows  easily  from  (5.6),  (5.7)  and  (4.16) 

that  pj(xj,0)  and  p^(xj,0)  satisfy 

a(pj(Xj,0))  -  b(Pj(xj,0))=Xj  for  -oo <Xj  £  a(l )  -  b(  1 )  ,  >j 


afPjfxj.O)) +b(pj(xlf  0))  =Xj  for  a(  1 )  -  b(l )  <Xj  <oo 


l  (5.11) 


a(p2(Xj,  0))  -  b(p2(xj,  0))  =  Xj  for  -oo<Xj£0 


a(o2(x^,  0))  +  b(c2(xj,  0))  =Xj  for  0<Xj<oo 


>  (5.12) 


Bearing  (3.9),  (3.10)  and  the  second  of  (3.6)  in  mind,  we  now  define 
functions  Uj(Xj,x2)  ,  u2(xj,x2)  for  all  (xj,x2>  by 

1/2 

ua  =  sgnx2Da{2b(oa)[a(p^)  +  b(0Q)  -X]]} 
a  =  1 , 2  ,  no  sum  . 

Note  that  Uj=u2  on  the  exterior  of  the  envelope  id.  From  (5.11),  (5.12), 
(5.  13)  one  sees  that 


r  (5.13) 


1 


Uj(Xj,  0+)  -  Ujfxj,  0-)*^ 


Vo 


=  0 


for  -oo<Xj  £  a(l)  -  b(l)  , 
for  a(l )  -  b(l  )<Xj  <oo  , 


where 


(5.  14) 


Since  b(p^)>0  ,  the  fact  that  Xj  lies  between  the  Xj -intercepts  of  r(p  ) 
assures  that  the  radical  in  (5.  13)  is  real. 
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u2(x1#  0+)  -  u2(x1#  0- ) 


^  0  for  -oo <Xj  -  0  , 


=  0  for  0<Xj<oo 


(5.  15) 


Thus  Uj  is  continuous  at  points  on  the  x^-axis  to  the  right  of  the  left 
x j -inte rcept  of  F(l)  but  discontinuous  across  the  crack  XjS  0  ,  =  0  , 


as  well  as  across  the 


;nt  O^XjS  a(l)  -  b(l)  of  the  -axi_s  ahead 


of  the  crack.  In  contrast,  u^  is  discontinuous  across  the  crack  but 
remains  continuous  at  all  points  of  the  Xj-axis  ahead  of  the  crack-tip. 
These  observations,  together  with  the  remarks  following  (5.7)  concerning 
the  loci  of  points  of  discontinuity  of  and  p^  ,  permit  the  determination 
of  the  maximal  domains  of  smoothness  for  Uj  and  .  Let  £j  consist 
of  the  set  of  all  points  in  the  domain  £  which  lie  neither  on  the  line 
segment  0<XjS  a(l)  -  b(l)  ,  x2=0  nor  on  ^iat  arc  envelope  =?  which 

joins  the  two  cusps  by  way  of  the  point  Xj  =  a(l)-b(l)  ,  x^=  0.  Similarly, 
let  f>2  be  the  set  of  those  points  in  £  which  do  not  lie  on  that  arc  of  s£ 
which  joins  che  two  cusps  via  the  origin.*  Then  u^  is  three  times  con¬ 
tinuously  differentiable  on  £'^  .  Moreover,  one  can  show  by  direct  cal¬ 
culation  based  on  (5.13),  (4.  16)  (and  the  fact  that  P  j  >  are  zeros  of 

F)  that  ua  satisfies  the  differential  equation  (2.  13)  on  .  In  the  process 

of  carrying  out  this  calculation,  one  finds  that 


u  =-sgnx2p&. 


P  a(p&)  +  b(p^)  -  x1  y  /2 


2b(pa) 


on  £ 


(5.  16) 


It  is  understood  that  the  cusps  themselves  do  not  belong  either  to  £.  or 
to  . 

^  v 


/ 

#  /  V*  *  •  . 
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l~x]  -a(p^)  +  b(pQt)  1 1 

uo?,  2"  po  2b(p^)  J 


on  £ 


(5. 17) 


and  hence  that 


1 7u  =  p  on  £ 

'  a  a  a 


(5.  18) 


In  view  of  (5.8),  (5.9)  and  (5.  18),  one  concludes  that  the  differential 
equation  (2.  13)  is  elliptic  at  the  solution  Uj  at  those  points  of  £j  which  lie 
outside  r ( 1 )  and  hyperbolic  in  that  part  of  £j  which  lies  inside  T(l).  Simi¬ 
larly,  (2.  13)  is  elliptic  at  the  solution  u^  at  those  points  of  which  lie 
outside  both  s£  and  P(l)  fu^  and  Uj  in  fact  coincide  there),  but  hyperbolic 
at  any  point  of  2L  which  lies  either  in  the  interior  of  P(l)  or  in  the  interior 


of  £  . 


1-2  .  _2 


One  has  p  j  (x^ ,  x^)  =  0  2(Xj  >  x^)  ide  si  and  as  r  =  'x^  +  x2  -  0  , 


their  common  value  tends  to  zero.  In  fact,  from  (4.  10),  (4.  16),  (5.6)  and 
(5.7)  it  follows  that 


p](x1,x2)  =P2(Xj,x2)  ~  2' 1/2  r' 1/2  as  r-oo  ,  (5.19) 


and  hence  from  (4.  10),  (5.  13)  and  (5.  19)  that 


1/2  o  _ 

u  j(x^,  x2)  =  u2(xj,  x2)  ~  (2r)  sin-^  as  r-*00  .  (5  .  20) 


Thus  Uj  and  u2  satisfy  the  condition  at  infinity  (2.  15)  posed  as  a  part  of 
the  small-scale  nonlinear  crack  problem. 

From  (5.  17)  with  a  =  1  ,  one  obtains 


/ 

■  / 
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U1  2^1'0±)  =  P  0) 


'xj-  a(p  jfxj,  0))  4  b(p  jCjEj,  0)) 


2b( 


M  1  i  #  0)  ) 


I 


/2 


(5.21) 


But  according  to  the  first  of  (5.  11),  the  numerator  in  the  bracket  in  (5.21) 
vanishes  for  Xj<a(l)-b(l)  ,  so  that 

ui  2(xj,0±)  =  0  *  -oo<x1<a(l)-b(l)  .  (5.22) 

Thus,  in  particular ,  u^  satisfies  the  free  surface  condition  (2.  14)  along  the 
crack*  -oo<Xj<0  ,  x^O. 

Neither  u^  nor  u^  alone  can  serve  as  a  solution  of  the  small-scale 
nonlinear  crack  problem,  since  neither  is  smooth  enough  throughout  *£  . 

In  the  next  section  we  shall  show  how  they  can  be  pieced  together  to  construct 
the  desired  solution. 

It  may  be  remarked  that  a  third  solution  u^  of  (2.  13)  differing  from 
Uj  and  u2  only  at  points  inside  ,  could  be  constructed  with  the  help  of 
that  zero  of  F  which,  for  (XpX^)  inside  JL  ,  lies  between  Dj  and  . 

We  shall  not  need  this  solution  in  what  follows,  and  so  we  do  not  investigate 
its  features  here. 

6 .  Solution  of  the  small-scale  nonlinear  crack  problem. 

Here  we  construct  a  continuous,  piecewise  smooth  solution  u  of  the 
boundary-value  problem  (2.  1 3 ) - ( 2 .  15);  u  will  coincide  with  Uj  on  a  subset 
of  H  and  with  u^  on  the  remainder  of  ?S  .  We  first  observe  that,  because 
of  the  discontinuities  in  u^  and  u^  ,  u  will  be  required  to  coincide  with  u^ 
at  those  points  of  2f  which  lie  near  that  branch  of  £  which  joins  the  cusps 

*  The  solution  also  has  vanishing  normal  derivative  at  x^=  0±  for 
-oo<x1<0. 
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by  way  of  the  origin,  while  u  and  u2  will  coincide  near  the  remaining 
branch  of  =£,  as  well  as  in  the  vicinity  of  the  segment  0  <Xj  <  a(  1 )  -  b(  1 ) 
of  the  3^-axis.  This  suggests  that  we  attempt  to  "match”  Uj  and  u2 
across  two  as  yet  undetermined  curves  «/+ and  J  (Fig.  8)  that  are  sym¬ 
metrically  situated  with  respect  to  the  x^-axis,  start  at  the  origin,  termi¬ 
nate  at  the  cusps,  and  otherwise  lie  wholly  within  the  interior  of  £.  .  This 
matching  must  assure  the  continuity  of  u  and  of  the  traction  across  the 
two  curves  to  be  found.  We  shall  show  that  there  are  unique  shocks  ® 
and  with  the  desired  properties.  As  either  of  these  curves  is  traversed, 
u  and  its  tangential  derivative  remain  continuous,  while  the  normal  deriva¬ 
tive  of  u  jumps  -  despite  the  prevailing  continuity  of  traction  across  J 
and  y  . 

Momentarily  taking  for  granted  the  existence  of  and  J  ,  let 
stand  for  the  set  of  those  points  in  7  which  lies  strictly  outside  the 
closed  curve  formed  by  J  ,  y  and  that  branch  of  which  connects  the 
cusps  byway  of  the  point  Xj=a(l)-b(l)  ,  x,=  0  (Fig.S),  while  *■  consists 
of  the  interior  of  this  closed  curve  (  f-  is  shaded  in  rig.  81.  Define  u  on 
7  by  first  setting 

*(1) 

on  £  , 

(6.  1) 

on  *<2)  , 

and  then  extending  the  definition  of  u  to  all  points  of  J  ,  J  (except  the 
origin)  and  the  right  branch  of  by  continuity.  It  is  clear  from  (5.20) 

that  u  satisfies  the  condition  (2.  15)  at  infinity  and  from  (5.22)  that  it 


i 
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also  satisfies  the  boundary  condition  (2.  14)  on  the  crack  faces.  Finally, 
u  is  smooth  and  satisfies  the  differential  equation  (2.  13)  on  £  except  at 
points  of  a^  +  and  J We  shall  ultimately  show  that  u  remains  bounded 
as  7-  0  . 

It  therefore  remains  only  to  establish  the  existence  of  curves 
and  ^"withthe  required  properties.  We  consider  «^+only;  the  existence 
of  a^'then  follows  from  symmetry. 

Suppose  for  the  moment  that  there  is  such  a  curve  and  let 

(x.,x2)  be  a  point  on  it.  Since  lies  on  or  inside  =£  ,  one  has  from  (5.  10) 
that 

0<o1(x1,x2)^  Pc  ,  p2(xlfx2)apc  ,  (xj,x2)6»;  (6.2) 

Since  u  as  given  by  (6.  1),  (5.  13)  is  to  be  continuous  at  (x^,x2)  on  J  , 
it  follows  that 


o2[2b2(a2+b2- Xj)]  1^2=pj[2b1(a1+b1- Xj)] 1//2  on  ,  (6.3) 

where  we  have  set 

bQ=b(oa)  ,  aa=a(Pa)  »  a=l»2  .  (6.4) 

Since  Oj  ,  c2  are  zeros  of  F(p;x^,x2),  one  has  from  (4.16)  that 

(X)- al>2+ x2  =  b2  ,  (xra2)2+x22  =  b2  on  J  (6.5) 

If  we  define  g  by 

g(p )  =  p2b(p )[ a(p)  +  b(p )]  ,  p > 0  ,  (6.6) 
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the  result  of  squaring  both  sides  of  (6.3)  and  solving  for  x^  may  be  ex¬ 
pressed  as 


-  g2~gl 
xl"  ,  2  2 

b2°2  "  blPl 


on  J  , 


whe  re 


ga=g(P0)  ,  «  =  1.2 


Using  (6.  7)  in  the  first  of  (6.  5)  yields 


(6.7) 


(6.8) 


on  .  (6.9) 

Substitution  from  (6.7),  (6.9)  into  the  second  of  (6.5)  then  leads  after  some 
algebra  to  the  conclusion  that  the  values  of  pj  and  p.,  on  the  shock  must 
be  related  by 

$(p1,D2)  =  0  on  ,  (6.10) 

where  i  is  defined  for  all  Oj>0  ,  p2>0  by 

2  2 

«  (p  P2)  =b1P  bj- a2+ b2)  -  b2P2(a2+ b2- aA+ bj)  ,  (6.11) 

(recall  (6.4)). 

The  issue  of  establishing  the  existence  of  the  solution  (6.  1)  of  the 

^Since  we  are  interested  only  in  p  j  >0  ,  p2>0  ,  (6.  7)  is  equivalent  to  (6.3) 
and  leads  to  no  extraneous  results.1 


-2  2 

*2=br 


go-  g- 


b2C  2  *  blpl 


2  '  al 


-37- 


small-scale  nonlinear  crack  problem  is  a  matter  of  the  reversibility  of 
the  argument  leading  to  (6.  10),  (6.  11).  We  shall  establish  below  the 
following  six  claims: 

(i)  for  every  o.  €(p^,p  ]  ,  there  is  a  uniaue  root  p,=  c,(p.)^p 

iuc  *  c  c-  1  c 

of  (6.  10); 

(ii)  setting  P ^  =  7 g(P  j)  in  (6.7),  (6.9)  supplies  a  parametric  rep¬ 
resentation  of  a  curve  •S+  with  tracing  parameter  p.; 

J.  *" 

(iii)  as  p.-*Pq+  ,  the  corresponding  point  on  J  approaches  the 

origin,  while  as  »  the  associated  point  on  J  moves 

to  the  upper  cusp; 

(iv)  except  for  its  endpoints,  lies  in  the  interior  of  =£  ; 

(v)  (6.3)  holds  on  »,+  (except  at  the  origin),  so  that  u  as  defined 
by  (6.  1)  is  continuous  except  at  the  origin; 

(vi)  the  traction  is  continuous  across  (except  at  the  origin). 

In  order  to  analyze  the  shock  equation  (6.  10),  we  require  a  further 
property  of  b(p)  as  well  as  some  results  pertaining  to  the  function  g(p) 
defined  in  (6.6).  From  (4.1),  (4.2)  and  (4.4)  one  finds  that 

[p2b(p)]'=  -p2a'(p)>0  ,  p  >0  ,  (6.12) 

2 

so  that  p  b(p)  is  an  increasing  function  of  p  .  Since  by  (4.  10), 

2 

p  b(p)-l/2  as  p-0+  ,  it  follows  that 

p2b(p)>-|  ,  P>0  .  (6.  13) 

The  asymptotic  behavior  of  g(p)  can  be  determined  from  (4.  10), 
(4.6)  or  (4.7),  and  (6.6)  as 

g(p) - 4 

4p 


as  p  -•  0+  , 


(6.  14) 


I 
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l 

1 


and 

2 

g(p) - as  p-  +  oo  if  T(X)>0  ,  (6.15) 

2t 

ao 


while 

g(p)  ~ -^y-p2^  '  €  ’  as  p-  +  oo  if  1^=0  .  (6.16) 

2h  c 

The  material  constants  t  ,  h  and  e  first  appeared  in  (1.  16);  note  that 

CD 

0<e  <  1  . 

From  (6.  6),  (6.  12)  and  the  first  of  (4.  19)  it  further  follows  that 


g'<P)  =  -P 


V(p)e  j(d) 


>0 


(6. 17) 


According  to  (6.6),  (6.14),  (6.17),  (4. 4)  and  the  graph  of  e^p)  in  Fig.  6, 

g(p)  must  decrease  from  too  at  0  =0+  to  a  positive  minimum  at  that 

value  of  p  at  which  e^p)  vanishes,  thereafter  increasing  steadily  with 

p  .  According  to  (6.  15),  (6.  16)  g  rises  to  a  finite  limit  as  p  -too  if 

t  >0  ,  but  tends  to  +oo  as  p-*+ao  if  t  =0  . 
ao  oo 

If  t  >0  ,  there  is  a  unique  number  p  >0  such  that 

CD  '  CD 

2 

g(poo)  =  S(+0°)  =  •  (6.18) 

2t 

CD 

With  the  help  of  a  cclculation  too  long  to  be  included  here,  it  is  possible  to 
show  that 


P  <p„ 

OO  ^  0 


(6. 19) 
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Pg  being  the  zero  of  a(p)-b(p)  (see  (4.15)).  Thus 


g(p)<g(oo)  =  for  ^  Too>0 

2t 

CO 


A  sketch  of  the  graph  of  g (p )  is  shown  in  Fig.  9. 

We  now  turn  to  the  function  $(P,»P2)  introduced  in  (6 .  1 1 ), 
which  it  is  immediately  noted  that 

«(P1,PJ)=0,  $(P2,  Pj)  =  (P  1#  p2)  . 


We  are  interested  in  the  behavior  of  5(Pj,02)  as  a  function  of  p0 
0<p2<ao  ,  when  is  held  fixed  at  a  value  not  less  than  p^.  A 

forward  calculation  based  on  (6.  11),  (6.  12)  and  the  first  of  (4.  19) 
the  identity 


b2fr(prD2)  =  a2[Y(prp2)'?(prp2)]  ’  Di>0  -  p2 


>0 


with  Y  given  by 

T(P1,P2)  =  (b2P2- b1Pj)e1(p2)  +  g1-g2  ,  Pj>0  ,  P2>0  , 

where  ga=g(PQ),  g  is  defined  by  (6 . 6 ),  e^p)  is  given  in  (4.  19), 
have  set  a'2=a'(p2).  From  (6.23),  (6.22),  and  (6.21)  it  follows  in 
that 


Y(pj,p1)=0  ,  Pj>0 


(6.  20) 

from 

(6.21) 

t 

straight- 

furnishes 

(6.22) 


(6. 23) 

and  we 
particular 

(6.24) 


and 
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(p  . ,  p  , )  =  0 

1  1 


Pj>0  . 


(6.25) 


Differentiating  (6.23),  (6. 22)  with  respect  to  c>2  ,  setting  P2=pj  ,  and 
making  use  of  (6.  12),  (6.  17),  (4.20)  yields 

a2* 

— |(P1,P1)=0  ,  p  j  >  0  .  (6.26) 

9e2 


Finally,  a  further  differentiation  gives 


lit..  .  'Pi^P 

3  (P  i ,  P  i )  -  b  Cj 

ap2  i 


(6.27) 


where  a'  =  a'(o  j )  ,  c^=c(Pj)  ,  and  the  function  c  is  defined  in  (4.21).  In 
view  of  the  one-cusp  condition  (4.24),  (6.27)  and  (4.4)  furnish 


83t 


ap 


3  *p1'd1H 


<0  , 

o  <  p , <p^ 

I  c 

=  0  , 

Pl  =  Pc 

(6.28) 

O 

A 

P>PC  . 

Making  use  of  (6.21),  (6.25),  (6.26),  (6.28)  and  the  appropriate  Taylor 
approximation  for  §(Pj,p2>  ,  one  can  thus  determine  the  local  behavior 
of  f(p  j,P2)  as  a  function  of  p2  near  the  zero  P2=pj  5 


$(Pj,P2)  decreases  with  p2  at  P2=Pj  if  0<p^<| 

t(Pj,P2)  increases  with  p2  at  P2=Pj  if  pi>pc 
From  (6.23),  one  finds  that 


>  (6.29) 
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9y  (p  p  )  2  2 

a*  :|Vi-tilil‘zcz  •  pi>0-p2>0  - 


(6.30) 


where  c,=  c(p7)  and  c(p)  is  given  by  (4. 21).  Invoking  the  one -cusp  condition 

(4.24),  one  sees  that  the  graph  of  y  (p  j ,  p  ^)  vs  p2  bas  horizontal  tangents 

only  at  p,=  p,  and  p^=p  .  Differentiation  of  (6.  30)  with  respect  to  p  ? 

7  21  2  c  ^ 

gives,  at  p  2=  P  j  » 


S2y  .  2.  ,  .2  n 

— jlp  p  l  =  -p  (a^  Cj  .  P  j >0  , 

3p2 


(6.31) 


so  that,  by  (4.  24), 


>0  ,  o<P 1<pc 


7T(pi'piM  =  0  ’  prpc 


<0  ,  Pl>pc  . 


(6.32) 


Next  we  consider  the  asymptotic  behavior  of  $(p^,p2)  and  Y(Pj,p2) 
as  p  2  +oo  and  as  C  2  ®4’  ^or  P  i  ^  P  q  •  ^  follows  easily  from 

(4.6)  or  (4.  7)  and  (6.  11),  (6.23),  (4.19),  (6.6),  (6.  18)  that  as  P2-°°, 
one  has 


for  t  >0: 
oo 


*  (P  1 '  P  2>  ~  ~2a  ^al "  bi)p  2>  0  if  pl>p0  ' 


(6.33) 


$(pQ,P2)  ~  g(pQ)  -  g(pQO)<o 


(6.34) 


7(p  ,,P,)  ~  g(p  ])  -  g(Poo><0  if  P  P n  ' 


1  K0 


(6.35) 


*  \  \ 

+  /  -&•  £•  "■  « 

/  •• 
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fo  r  t  =  0  : 
oo 


$(pi'p2)~‘^(arbi)p2’€>0 


if  p  1>P0  • 


(6.36) 


*(p0'p2}~*  T7rp2(1’€)<0 
£ 


(6.37) 


7(p1,p2) 


— - p^(1*s)<0  if  p,>0  . 

2h(2  -  e )  C 


(6.38) 


On  the  other  hand,  the  small-p^  behavior  of  $  and  Y  can  be  deduced 
with  the  help  of  (4.  10):  as  P2“*0  +  , 


5  (p  j »  P  2) 


(bi'i 


Pj>0 


(6.39) 


y(p  i'p2) 


2dp 


4<bl 


p?-4)>o 


Px>0 


(6.40) 


We  shall  now  make  use  of  all  of  the  foregoing  properties  of  J  and 
Y  in  order  to  describe  the  graphs  of  Y(pj,p2>  and  $(Pj,P2)  as  functions 
of  p2  ,  0<p2<oo  ,  with  pj  fixed,  pQ£  p  pc  .  Suppose  first  P0<P1<PC  * 
From  (6.40)  one  sees  that  Y(pj,p2)  decreases  from  +oo  at  p2=0+  until 
it  reaches  its  zero  at  p2=p^  (see  (6.24)),  at  which  point  it  has  a  horizontal 
tangent,  according  to  (6.30).  As  p2  increases  from  Pj  ,  Y(pj,p2) 
increases,  reaching  a  maximum  at  p  =pc  ,  thereafter  declining  steadily. 

If  t  >0  ,  Y (p. , ,  p  , )  approaches  a  finite  negative  limit  as  p?-  +oo  ,  by 
(6.35).  If  Too=0  ,  (6.38)  shows  that  Y(p  ^ ,  p  2) -• -cd  as  p  2  -*  +oo .  The 
graph  of  Y(pj,p2)vs.  p2  for  Po<pi<pc  *s  shown  as  the  dashed  curve 


/ 


S&' 


Vi  , 
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in  Fig. 10. 

According  to  (6.22),  the  graph  of  i(pltp  )vs.  p2  can  have  a 
horizontal  tangent  only  where  it  crosses  the  graph  of  Y(pj,p2)  vs.  p^.  By 
(6.39),  (6.22),  ?(fj,p2)  decreases  from  +co  at  p2=0+  ,  remaining  less 
than  Y(pj,p2)  until  it  vanishes,  together  with  Y(Pj,p2>,  at  P  2=  p  1  ' 
where  there  is  a  horizontal  tangent  (see  (6.25)).  By  the  first  of  (6.29), 

5  (p  i »  P  2 )  *s  negative  and,  by  (6.22),  continues  to  decrease  as  p2  increases 
from  pj  until  its  graph  crosses  that  of  Y(pj,p2>  with  a  horizontal  tangent. 
As  p^  continues  to  increase,  }(pj,p2)  then  increases  monotonically  to 
+co  (see  (6.33)  or  (6.36)),  crossing  the  p^-axis  once  and  only  once  at 

P2=P2^V*  graph  of  f  (p  ^ ,  p  2>  vs .  p  2  for  p0<pl<pc  is  shown  as 

the  solid  curve  in  Fig.  10. 

For  Pj  =  Pq  or  p  ^  =  p  c  ,  it  is  easily  shown  that  $(p  ,p7)  vanishes 

if  and  only  if  p,  =  p,  ,  whether  t  >0  or  r  =0  . 

c  1  oo  co 

Thus  for  Pq<pi^pc  >  there  is  a  unique  root  p2=p"2(pj)2p^  — 

(6.  10),  whether  Too'>0  or  Tao=®  •  an£l  claim  (i)  is  established.  A  detailed 
study  of  the  behavior  of  p2(p  j)  shows  that  p"2(pj)  decreases  monotonically 
as  Pj  increases,  and  that 


Pz(pl)-+°o  as  p  j  ~*P  q+  • 
P^P  i)_,pc+  as  pl"pc*  • 


r  (6.4i) 


In  fact,  the  asymptotic  behavior  of  p^fpj)  as  pl“,po+  can  ^  shown  to  be 
given  by 
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1 


P2(pl)~  P0 


■  2,  .  2' 

T  (  C  «  )  -  T 

ao 


Ta>T'<l’ol 


(Pl-P0) 


-1 


if 


T  >0 
ao 


(6.42) 


Here  pg<l  *s  th®  unique  zero  of  a(p)-b(p)  (see  (4.15)),  T(k)  is  the 

shear  stress ^  in  simple  shear  at  an  amount  of  shear  k  ,  and  t  is  the 

oo 

ultimate  stress  in  simple  shear. 

If  Too=0  >  (6.42)  must  be  replaced  by 


P2(P 


r  2-  v 

poT  (po} 

€hT/(p0) 


1 

e 


(Pi-Pq) 


1 

€ 


(6.43) 


here  h  ,  e  are  the  material  constants  introduced  in  (1.  16),  (1.  17). 

If  one  sets  p^=p'-)(p  j)  in  (6.7),  then  Xj  is  given  as  a  function  of 

c.  for  p  .  <  p  .  ^  p  .  In  order  to  show  that  the  same  substitution  in  (6.9) 
10  1c 

furnishes  as  a  function  of  p  ^  ,  it  is  first  necessary  to  establish  that 

the  right  side  of  (6. 9)  —  with  p  2=  P  2(p  j )  ~  is  positive .  Let 


If 


p  T  P  2 


^(Pj»P2)~bj 


2  I  g2‘gl 


2  2  *  a  1 
<b2p2  ~  blp  1 


satisfy  (6.  10)  with  $  given  by  (6.  11), 


2>0  ,  p  ^ > 0  .  (6.44) 

it  is  possible  to  show  that 


b  b  22 

z(p  i ,  p  2>  =  1 1  1  %  2(ar  a2+brb2)(ar  a2'bi+b2) 

b2p  2  ’  blP  1 


(6.45) 


When  P  2>P  j  »  one  has  W  (4.  4)  that  a^  -  a2+  bj  -  b2>0  ,  so  that  the  sign  of 


*Note  that  T(k)  is  the  dimensional  shear  stress  (as  distinguished  from  T(k) 
which  is  nondimensional). 
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Z(p  i .  p  2)  is  the  same  as  that  of  a^- a^- b^+ ,  provided  (6.  10)  holds .  To 
determine  the  latter  sign,  we  first  use  (6.11),  (6.6),  and  (4.19)  in  (6.23) 
to  show  that,  if  p  ^ ,  p  2  satisfy  (6.10),  then 


v(p  l*  p2)  "  a ^  (a2*  b2,(b2p2  "  blp  1 )  *  b2p  2(al '  a2‘  bl+b2) 


(6.46) 


But  according  to  Fig.  10,  y(p  p  ^(p  ^))<0  .  Moreover  a^cO  by  (4.4), 

aa‘b2<0  if  P2  =  ^VP1^Pc>1  hy  {4’  5)’  311(1  b2p2"blPl>0  if  P2>P1 
by  (6.  12).  When  these  facts  are  used  in  (6.46),  one  concludes  that 

a  j  -  a,  -  bj  +  b^  >0  when  P2=P^P1^'  and  hence  that  Z(pj,p '^(p^))>0, 

Pq<Pj^Pc>  It  then  follows  that  the  right  side  of  (6.  9)  is  positive  when  p  2=  p"^(p  j )> 

and  hence  (6.  7),  (6.  9)  with  p2=p"2(p  j),  p  Q<  p  p  c  supply  parametric  equations 

of  a  curve  for  x->>  0  (or  J~  if  xo<0),  thus  e stablishing  claim  (ii). 

U  L. 

Equation  (6.41)  shows  that  P*2(P  j ^  +oo  as  Pj-*Pq+  ,  and  thus  by 
(4.  6)  or  (4.  7),  (6.  15)  or  (6.  16),  and  (6.  7),  it  follows  that  Xj  -0+  on  J* 
as  p  j  -»Pq+  .  A  similar  argument,  making  use  of  (4.  15)  as  well,  shows 
that  X2-*0+  on  t/+  as  Pj-»Pq+. 

On  the  othe r  hand,  as  p  j  ~ ’  P  c ~  »  (6.41)  asse rts  that  P  2^p  l )**P  c+  • 

(6.17),  (6.12),  (4.19),  (6.7)  and  L'Hospital's  rule  then  show  that 

Xj-*ej(pc)  on  .  It  then  follows  from  (6.9/  and  (4.  19)  that  -'Je^p 

on  as  pj-»p  -.  Thus  (x ^ in  approaches  the  upper  cusp 

of  sd  as  P  ,  -*p  -  ,  and  claim  (iii)  is  confirmed. 

1  c 

Suppose  that  (x^x^)  is  an  interior  point  of  a/+  corresponding  to 
p  j  ,  pQ<p  j<Pc  •  Then  (6.  10),  (6.  11),  (6.  7),  (6. 9)  and  (4.  16)  imply  that 
Pj  and  ^(Pj)  are  both  zeros  of  Ffpjx^x^)  ,  so  that  (x^,^)  is  a  point 
at  which  F  has  more  then  one  zero.  By  Table  1  of  the  preceding  section, 


'  /  • 
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(x  ,x2)  cannot  lie  outside  at.  If  (Xj,x2 )  were  on  at,  then  pj  and 
p"2(p  would  both  be  double  zeros  of  F  (see  (4.  17))  ,  which  is  easily 
shown  not  to  be  the  case.  Thus  (xj,x2)  lies  inside  at ,  as  asserted  in 
claim  (iv). 

It  can  be  shown  directly  that  (6.  10),  (6.11),  and  (6 . 7)  —  with 
p  2=  P"2(P  i )  »  Po<pl^pc  _  (6*  3),  and  hence  by  (6.  1 ),  (5.13),  u  is 

continuous  at  all  points  of  except  the  origin.  This  establishes  claim  (v). 

To  verify  that  the  traction  is  continuous  across  the  shock,  it  is  nec¬ 
essary  first  to  determine  the  slope  of  .  If  one  sets  p  ^ =  P*2 (P  j )  *n  (6.10) 
and  then  differentiates  with  respect  to  p  one  can  determine  p"2(p  j)  in 
terms  of  ,  p^p^P})  with  the  help  of  the  appropriate  properties  of  a(p  )  , 
b(p )  and  i(p  j,P2)  •  To  find  the  slope  of  a/+  ,  it  is  then  only  necessary  to 
differentiate  (6.7)  and  (6.9)  (with  p2=p"2(p  j))  with  respect  to  Pj  and  sub¬ 
stitute  for  p" ^ (p  j )  the  value  determined  as  outlined  above.  The  resulting 
formula  for  the  slope  can  be  simplified  by  using  (6.  10),  (6.  11)  to  eliminate 
p^,p2  where  they  appear  explicitly.  The  final  form  for  the  slope  is  then 
found  to  be 

’(aryyb2)(«2- vyM  1/2 

.(wyty(v*2+y  b2>. 

on  ,  p2=p'2(p1),  P0<P1<PC  •  (6.47) 

We  note  in  passing  that  the  monotonicity  properties  of  a(p  )  ,  b(p )  and 

2 

p  b(p)  established  in  (4.4),  (6.  12)  permit  one  to  conclude  that  the  slope 


_  S  _  Vb2 

m~  dXj  '  bj+b2 
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of  J*  is  positive  at  each  of  its  interior  points.  1 
The  vector  N  with  components 


N^=  -m(  1  +  m  ) 


,  N^=  ( 1  +  m  ) 


(6.48) 


in  which  m  given  by  (6.47),  is  a  unit  normal  vector  of  J  .  By  (2.6), 
(2.  7),  the  only  nonvanishing  component  of  traction  on  a  curve  in  l  with 
normal  N  is  the  x^  -component,  which  is  given  by 


t  =  T3pNp=2W'(3+  |Vu|2)u>(3Np 


(6.49) 


By  (1.6),  (3.2)  and  (4.1),  this  may  be  written  as 


—f— G  6N6  '  P  =  'V-51 

2p  d(p  )  'P  P 


(6.50) 


The  components  u  .  ,  u  ?  of  7  u  are  discontinuous  across  ,  being 
given  by  (5.  16),  (5.  17)  and  (6.  1)  on  the  side  of  associated  with  p^  , 
a  =  1,2  .  The  limiting  values  of  u  .  ,  u  -  on  the  p  -side  of  can  be 

,1  p  L.  Of 

expressed  entirely  in  terms  of  p  ^  by  substituting  for  x^  in  (5.  16), 

(5.17)  from  (6.7).  Since  Nj,  are  expressed  solely  in  terms  of  Pj.P^ 
by  (6.48),  (6.47),  the  limiting  tractions  tj  ,  t^  on  the  two  sides  of 
can  then  be  found  in  terms  of  p^,p ^  alone  by  (6.50).  In  the  resulting 
formulas  for  tj  and  t^  ,  aPPear  explicitly  as  well  as  in  the  argu¬ 

ments  of  a  and  b.  It  is  possible  to  use  the  shock  condition  (6.  10),  (6.  11) 
to  eliminate  the  explicit  appearance  of  p^,p ^  .  When  this  is  done,  one 

*Note  that  one  also  requires  the  fact  that  aj-a^-b^  b^>0  as  established 
in  the  discussion  following  (6.46). 
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finds  that  tj=t2  on  **+  »  confirming  the  continuity  of  traction  across 
and  thus  establishing  claim  (vi). 

Geometrically,  the  shocks  and  J  are  the  loci  of  points  of 
intersection  of  the  circles  T(p^)  and  T(p,)  as  Pj  varies  from  p^  to 
p  ,  p_  from  +oo  to  p  in  accordance  with  the  relation  p_=p,(p.)  . 

Cl.  C  L.  l.  I 

These  circles  are  shown  schematically  in  Fig.  8. 


7.  The  field  near  the  crack-tip.  Discussion 

The  local  behavior  near  the  crack-tip  of  the  shock  and  of  the 
displacement  and  stresses  can  be  determined  from  the  global  solution  to  the 
small-scale  nonlinear  crack  problem  constructed  in  the  preceding  section. 
We  consider  first  the  case  of  a  material  for  which  the  ultimate  shear  stress 
in  simple  shear  is  positive:  t  >0  •  To  determine  the  shape  of  near 
the  origin,  it  is  easiest  to  proceed  from  the  expression  (6.47)  for  the  slope 
m  .  Since  Pj-*Pq+  and  p^-^oo  on  »'+  as  the  origin  is  approached,  one 
can  determine  the  limiting  value  of  dj^/dx^  as  Xj-»0+  from  (6.47)  by 
using  the  asymptotic  results  (6.42)  and  (4.6),  along  with  (4.15),  (4.  1), 

(4.3)  and  (3.2).  This  leads  to 


lim 


T 

GO 


>0  , 


(7.  1) 


where  tq=t(Pq)  •  and,  by  (4.15),  (4.1)  and  (3.2),  p^  is  the  unique  number 
in  (0,  1)  such  that 


— - 

■■  >  is 

•J 

t  V  *  /!.  .  ,r 
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PnT(Pn) 


(7.2)' 


PnP  T(P)  P°  H0 


From  (7.  1)  it  follows  that  J  is  asymptotically  straight  at  the  origin: 


x-,  - - — - x,  on  »,+  as  x,  -0  ,  t  >0  . 

2  r 2 - T  1  1  00 

0  oo 


(7.  3) 


Thus  J  and  J  subtend  a  limiting  angle  29^  at  the  origin,  where 


a  .2  2-1/2 

tan  G  _  =  t  (t  -  t  ) 

0  oo'  0  ao 


To  find  the  limiting  value  of  ufxj.x^)  as  (Xj.x^)  approaches  the 
origin  from  within  the  domain  (i.e.  from  the  elliptic  side),  one  uses 

(6 .  1 )  and  (5 .  13 )  with  a  =  1  and  lets  p  ^  — •  p  ^  +  .  Making  use  of  (4.  15),  (4.1), 
(3.2)  and  (7.3),  one  obtains 


lim  u(x. ,  x_)  =  ±—  as  r-0,  G  _.  <  t  9  '  ^  tt  ,  r  >  0  , 
r-0  1  2  70  0  °° 


(7.  4) 


where  the  ambiguous  sign  is  chosen  to  be  the  sign  of  0.  Note  that  (7.4) 
confirms  that  u  is  discontinuous  at  the  origin. 

In  order  to  find  the  limiting  value  of  u  at  the  origin  from  the  hy¬ 
perbolic  domain  ,  it  is  necessary  first  to  observe  that  the  second  of 

(6.  5),  together  with  (4.  6),  implies  that 


P2(Vx2)~?i“£¥1  as  r-°'  Too>0  ' 

oo 


(7.5) 


One  can  use  (4.  15),  (6. 6),  the  second  of  (4.  1),  and  (6.  20)  with  p  =  p( 

to  show  that  t.>t  , 

0  ao 


'  /  ^  h  - 

I  '■  ■:*  •/. 
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Appealing  to  (6.  1),  (5.  13)  with  a  =  2  ,  letting  r-0  and  invoking  (7 .  5), 


one  is  led  to 


lim  u(x, ,  x0)  =  —  sin  9  ,  -6^2  0  2  ,  Too>0 


‘r  2'  t 


(7.  6) 


A  graph  of  the  limiting  value  of  u  vs.  9  ,  for  T00'>®  *  is  shown  in  Fig.  11. 
Note  that  (7.4),  (7.6)  together  show  that  u  is  bounded  near  the  crack-tip, 
as  originally  required. 

The  corresponding  limiting  values  of  the  stresses  can  be  found 

from  (1.6),  (1.16),  (2.6),  (5.  l6)-(5.  18)  in  an  analogous  way.  The  cylindrical 

components  of  shearing  stress  t  =  T  3  jCos  9  +  6  and 

t  .  =  -r  ,  sin  9  +  T. -cos  9  turn  out  to  satisfy 
9z  31  32 


lim  t 


o  ,  -eo<9<0o 


±TqCOS  0  ,  9q<  I  0  !  ^  TT  , 


T  >0 

OO 


(7.7) 


lim  t, 


02  9  <  9, 


Tq | sin  9 


T  >0 

CO 


9q<  !  6  ;  ^  tt  , 


(7.8) 


The  ambiguous  sign  in  (7.7)  is  opposite  to  the  sign  of  9.  Graphs  of  these 
limiting  stresses  (t^>0)  are  given  in  Fig.  12,  in  which  9q  was  chosen 
to  be  greater  than  tt/4. 


^A  more  refined  calculation  shows  that  T^^Ofr)  for  -0q<0<9q. 


L  <n'<  ■  •  v,  •, 

'  / 

i  •  . 
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The  axial  stress  t  is  given  asymptotically  by 


'u  c°s  8  e 

r  '  S0<9<90 


33  ^ 


as  r  -  0  ,  t  >0 

oo 


(7.9) 


T  D  0  <  ®  S  ■ 

0P0  '  -0  1  '  1 


The  limiting  behavior  as  r-*0  of  u  ,  t  ,  and  is  quali¬ 

tatively  similar  to  the  small--?  behavior  of  these  quantities  as  determined 
in  [  l]  for  the  small-scale  nonlinear  crack  problem  in  the  case  of  an  in¬ 
compressible  material  governed  by  the  shear  stress  response  (1.  13).  For 
the  displacement  and  the  shearing  stresses,  this  similarity  can  be  seen  by 
comparing  Figs.  11,  12  of  the  present  paper  with  Figs.  7,  8,  9  of  [  1  ] . 

Apart  from  the  presence  of  the  discontinuity  in  stresses  associated 
with  the  shock,  the  above  limiting  values  of  displacement  and  stress  also 
bear  a  qualitative  resemblance  to  corresponding  results  obtained  in  [2]  for 
the  "softest"  of  the  elliptic  power-law  materials  (the  case  n  =  1  /2  in  ( 1 .  1 1)). 

It  is  possible  to  confirm  directly  the  continuity  of  displacement  and 
traction  across  J  near  the  origin  by  making  use  of  (7.  3),  (7.  4)  and  (7.  6)- 
(7.8). 


For  the  case  of  a  material  with  t  =0  ,  similar  —  but  more  in- 

oo 

volved  —  arguments  which  we  shall  not  give  here  lead  to  the  appropriate 
crack-tip  approximations.  In  this  instance,  the  shocks  and  J  are 
tangent  to  the  Xj-axis  at  the  origin  according  to  the  formula 


-Cl!N1/V/,3c 


1/e 


\u 


1 


as  x  -•  0  + 


on 


t  =0 
co 


(7.  10) 
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Here  TQ=T(p0*  '  Pq  determined  by  (7.2)  ,  and  h  ,  e  are  the  material 

constants  governing  the  shear  stress  response  at  large  shear  (see  (1.  16)); 
recall  that  h>0  ,  0<e<  1  . 

The  limiting  value  of  the  displacement  ufxjjX^)  at  the  origin  from 

within  "2^  (i.e.  from  the  elliptic  side)  is  again  given  by  (7.4),  except 

that  now  9q  =  0  ,  since  and  J  are  tangent  to  the  Xj-axis  at  the  origin. 

On  the  other  hand,  if  the  origin  is  approached  along  a  curve  which  lies  in 

~(2) 

the  hyperbolic  domain  £  ,  (7.6)  must  be  replaced  by 


-l/2/uN1/e,l-l/s'2-e  nir^A1/e'I/2  sine  _-(l-e)/e 

u~e  t  2  - cos9+R(9) - pry  r 

h  '  e  [  cos  9  R(6)]  ^ 


as  (x, ,  x9)  (0,  0)  inside  f,  ,  t  =0  , 

1  4-  CO 


(7.  11) 


whe  re 


r  (2  -  e  )2  2.  4(l-e)"1/2 

R(S)=  2 — cos  3 - ‘ — 2 — 

€  e 


(7. 12) 


With  reference  to  (7.  12),  one  must  note  that  R(6)  is  real  for 
|  co  s  G|  ^  2(1  -  e)^^2/(2  -  e)  ,  and  only  small  values  of  |  S  |  are  relevant  in 
(7.  11),  again  because  and  of  are  tangent  to  the  x^-axis  at  the  origin. 

It  can  be  shown  from  (7.  10),  (7.  11)  that  u  remains  bounded  as  the 
origin  is  approached  from  the  hyperbolic  domain  along  any  curve  that  re¬ 
mains  between  and  J  .  Since  according  to  (7.4)  such  boundedness 
clearly  prevails  for  any  approach  to  the  origin  from  the  elliptic  side,  one 
concludes  once  again  that  u  is  bounded  near  the  crack-tip. 

The  choice  T(X)=  0,h  =  u,e-l/2  in  (1.  16)  furnishes  precisely  the 
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large-k  behavior  of  the  material  described  by  (1.  13)  and  treated  in  [l]. 
Although  this  material  is  not  included  in  the  class  of  materials  studied  in 
the  present  paper,  principally  because  of  the  lack  of  smoothness  at  k  =  1 
in  its  shear-stress  response  curve  (Fig.  2),  it  is  nevertheless  true  that 
setting  c  =  1/2  ,  h=u  in  (7.  10),  (7.  11),  (7.  12)  reduces  (7.  10),  (7.  11) 
precisely  to  the  corresponding  results  in  [  1  ] .  *  This  reflects  the  fact 
that  crack-tip  behavior  within  the  hyperbolic  region  is  determined  entirely 
by  the  large-k  behavior  of  the  shear-stress  response  curve.  It  is  also  true 
that,  if  one  uses  T(p)  as  given  by  (1.  13)  in  (7.2)  to  determine  ,  sets 
~0=T(p0)=UPo  »  and  uses  this  in  (7.4),  the  result  coincides  with  the  limit¬ 
ing  value  of  u  from  the  elliptic  side  at  the  origin  as  determined  in  [  1  ] . 

Crack-tip  approximations  for  the  stresses  in  the  case  t^=0  re¬ 
main  the  same  as  for  the  case  r  >0  insofar  as  limiting  values  from  the 

QO 

elliptic  side  are  concerned.  Approximations  appropriate  to  an  approach  to 
the  crack-tip  from  the  hyperbolic  domain  for  t  =  0  have  been  found  but 
will  not  be  given  here.  They,  too,  reduce  to  corresponding  results  in  [  1  ] 
when  e  =  l/  2  ,  h=u. 

As  a  point  on  the  shock  approaches  the  cusp,  the  shock  "weakens" 
in  the  sense  that  discontinuities  in  the  displacement  gradients  tend  to  zero. 
At  the  cusp  itself,  one  can  determine  the  field  quantities  in  terms  of 
Pc“the  value  of  p  at  which  c(p  )  in  (4.21)  vanishes.  A  typical  formula 
is  that  for  the  displacement  at  the  cusp: 


^ee  (4.  1),  (4.2)  and  (4.9)  of  [  1  ] . 


•  «  i,x 

'  / 

I  '  * 
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u  =u 
cusp 


■P  T(  p  )'x~  - 

c  c  dp \p 


T  Y 


1/2 


P=P. 


(7.  13) 1 


The  largest  resultant  shearing  stress  (t^  t^)*^  occurs  at 

points  of  the  circle  T(l)  and  has  the  value  t  (see  (1.  18)).  The  axial 

m 

stress  t  3  becomes  infinite  as  the  crack-tip  is  approached  from  within 
the  hyperbolic  domain. 

The  small-scale  nonlinear  crack  problem  appears  to  be  substantially 
more  complicated  if  the  one-cusp  condition  (4.24)  fails  to  hold.  If  there  is 
more  than  one  cusp  on  each  of  the  branches  of  the  envelope  si  ,  it  can  be 
seen  that  there  must  be  at  least  three  on  each  branch.  The  situation 
sketched  schematically  in  Fig.  13  would  seem  to  be  a  possible  one  as  regards 
an  envelope  si  with  three  cusps  per  branch.  Among  the  possibilities  sug¬ 
gested  by  Fig.  13  is  one  in  which  the  shocks  </+  and  J  are  forked,  as 
indicated  by  the  dashed  lines. 

In  [3],  [4]  it  was  proposed  that  elastostatic  shocks  should  be  subject 
to  a  certain  "dissipation  condition"  analogous  to  the  requirement  of  increas¬ 
ing  particle  entropy  across  a  gas-dynamic  shock.  For  the  shocks  constructed 
in  [l]  in  connection  with  the  small-scale  nonlinear  crack  problem  for  the 
special  material  (1.  13),  it  was  verified  that  this  dissipation  condition  was 
fulfilled  when  the  amount  of  shear  prescribed  at  infinity  is  a  quasi- statically 
increasing  function  of  time,  corresponding  to  lo  ad  in  g  (as  distinguished  from 
unloading) .  We  have  not  investigated  the  corresponding  issue  in  the  more 
general  setting  of  the  present  paper. 

*  Re  call  that  (1.  15)  implies  d/dp  (t  (p  )/p  )<0  . 
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